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Abstract 

The fraction of Z° — * bb events in hadronic Z° decays has been measured by the OPAL ex- 
periment using the data collected at LEP between 1992 and 1995. The Z° — > bb decays were 
tagged using displaced secondary vertices, and high momentum electrons and muons. Systematic 
uncertainties were reduced by measuring the b-tagging efficiency using a double tagging technique. 
Efficiency correlations between opposite hemispheres of an event are small, and are well understood 
through comparisons between real and simulated data samples. A value of 

cr(e + e _ — > bb) 

i? b = -r-r 1 i n = 0-2178 ±0.0011 ±0.0013 

er(e+e — > hadrons) 

was obtained, where the first error is statistical and the second systematic. The uncertainty on R c , 
the fraction of Z° — > cc events in hadronic Z° decays, is not included in the errors. The dependence 
on Rc is 

ARb = -0.056-^ 



Rb R c 

where AR C is the deviation of R c from the value 0.172 predicted by the Standard Model. The 
result for R\, agrees with the value of 0.2155 ± 0.0003 predicted by the Standard Model. 
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1 Introduction 



The partial width for the decay Z° — > bb is of special interest in the Standard Model. Electroweak 
corrections involving the top quark affect the Z° — > bb partial width, T b ^, differently from the widths 
for lighter quarks. As a result, the fraction 

r b E _ r(z° -» bb) 

Thad r(Z° — > hadrons) 

depends on the top quark mass, m t0 p, but has negligible uncertainty from the unknown Higgs boson 
mass and the strong coupling constant a s . The fraction r^/r^ad is a l so sensitive to various extensions 
of the Standard Model involving new particles such as additional quarks and gauge bosons, or the 
virtual effects of new scalars and fermions such as those expected in supersymmetric models |IJ. 
The quantity measured in this analysis is the cross-section ratio 

_ <r(e + e - — ► bb) 
b a(e + e~ — > hadrons) 

at the Z° resonance. This differs from the partial width ratio T^/T^^ because of the additional 
contribution from photon-exchange diagrams. These have been evaluated within the Standard Model 
using the program ZFITTER 5.0 0| which predicts that i?b is 0.0002 smaller than r^/r^ad- By 
convention, Z° — > bb events where an additional qq pair is produced via gluon splitting are included 
in the numerators of the definitions of r^/Thad an d Rb- The small number of events where the only 
bb pair is produced via gluon splitting, rather than directly from Z° decay, is not included in the 
numerators. The interference between these two processes, and the effect on the measurement of R^, 
are expected to be negligible |J. 

This paper supersedes the previously published OPAL measurement § of the fraction of bb events 
in hadronic Z° decays. The measurement is improved by employing higher performance vertex tag- 
ging using the upgraded silicon detector with two coordinate readout installed in 1993 ||, higher 
performance electron identification, and by including the data taken in 1995. 

The paper is organised as follows. The analysis method, based on the double tagging technique, 
is described in the next section. The OPAL detector, the selected event sample and the Monte Carlo 
simulation are reviewed in Section ||. The b-tagging methods are discussed in Sections ||| and ||. The 
result is presented in Section ||, with systematic errors being evaluated in Sections [?] and ||. A summary 
is given in Section ||. 

2 Analysis Method 

The analysis method is based on the double tagging technique. Each selected hadronic Z° decay 
event is divided into two hemispheres by the plane perpendicular to the thrust axis and containing 
the interaction point. A b-tagging algorithm is then applied separately to each hemisphere, and the 
number of tagged hemispheres iVt and events where both hemispheres are tagged iVtt counted in the 
sample of iVhad selected hadronic events. These quantities are related by: 

N t = 2N had {e h R h +e c R c + e uds (1 - R h - R c )}, (1) 
Ntt = N had {C b (e b ) 2 R h + C c (e c ) 2 R c + C uds (e uds ) 2 (1 - R h - R c )}, (2) 

where e b , e° and e uds are the tagging efficiencies for hemispheres in bb, cc and light quark (uds) 
events, and C b , C° and C uds describe the tagging efficiency correlation between the two hemispheres 
in events of each flavour. For a useful b-tagging algorithm, e b is much larger than e c and e uds . The 
correlation C b is defined by C b = e bb /(e b ) 2 , where e bb is the efficiency to tag both hemispheres of a 
bb event. Deviations of C b from unity account for the fact that the tagging in the two hemispheres is 
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not completely independent, there being a small efficiency correlation between them for both physical 
and instrumental reasons. These correlations are also present in cc and light quark events; however 
their effect on the measurement is less than 10 because the double tagging efficiencies for these 
events are very small. The correlations C c and C uds are therefore set to unity. The values of and 
e b can then be obtained by solving equations |l] and ^ and only the values of e c , e uds and C h need to 
be input from Monte Carlo simulation. This technique avoids the need to input e° from simulation, 
which severely limits the precision of single tag measurements of i?b- 

In this analysis, two b-tags are used: a secondary vertex tag with very high b-purity, and a high 
momentum lepton tag with somewhat lower b-purity and efficiency. A hemisphere is taken to be tagged 
if it is tagged by either one or both of the secondary vertex and lepton tags. As the performance of the 
tagging algorithms changes significantly between the different years of data taking due to changes in 
the detector configuration, the equations are solved separately for e° and i?b for each year. The values 
of i?b are then combined, taking into account correlated and uncorrelated systematic uncertainties. 
The charm tagging efficiencies e c , light quark tagging efficiencies e uds and bb correlation coefficient C b 
for each year are input from Monte Carlo simulation, and are thus sources of systematic error. The 
relative sizes of these systematic errors and the statistical error depend on the purity and efficiency 
of the b-tag, and the tag cut is adjusted to minimise the total resulting error. The only other sources 
of systematic error are the value of R c (which is fixed to its Standard Model value) , and the hadronic 
event selection. 

3 Data Sample and Event Simulation 

The OPAL detector has been described in detail elsewhere 0,0]. Tracking of charged particles is 
performed by a central detector, consisting of a silicon microvertex detector, a vertex chamber, a 
jet chamber and z-chambers[]. The central detector is positioned inside a solenoid, which provides a 
uniform axial magnetic field of 0.435 T. The silicon microvertex detector consists of two layers of silicon 
strip detectors; the inner layer covers a polar angle range of | cos#| < 0.83 and the outer layer covers 
| cos#| < 0.77. This detector provided <p coordinate information in 1992, and was upgraded to provide 
both 4> and z coordinate information from 1993. In 1995, a new detector geometry was installed with 
smaller gaps between the wafers in 0, increasing the fraction of tracks with silicon hits from both 
layers. The vertex chamber is a precision drift chamber which covers the range | cos#| < 0.95. The jet 
chamber is a large-volume drift chamber, 4.0 m long and 3.7 m in diameter, providing both tracking 
and ionisation energy loss (d-E/dx) information. The z-chambers provide a precise measurement of 
the z-coordinate of tracks as they leave the jet chamber in the range | cos 6*| < 0.72. The coil is 
surrounded by a time-of-flight counter array and a barrel lead-glass electromagnetic calorimeter with 
a presampler. Including the endcap electromagnetic calorimeters, the lead- glass blocks cover the range 
| cos#| < 0.98. The magnet return yoke is instrumented with streamer tubes and serves as a hadron 
calorimeter. Outside the hadron calorimeter are muon chambers, which cover 93% of the full solid 
angle. 

The data used for this analysis were collected from e + e~ collisions at LEP during 1992-1995, with 
centre-of-mass energies at and around the peak of the Z° resonance. Hadronic events were selected by 
the algorithm used in giving a hadronic Z° selection efficiency of (98.1±0.5) % and a background of 
less than 0.1%. The thrust value T and thrust axis polar angle Ot were then calculated using charged 
tracks and electromagnetic clusters not associated to any tracks. To ensure a well defined thrust axis 
direction within the acceptance of the silicon microvertex detector, the thrust value and axis direction 
were required to satisfy T > 0.8 and | cos#t| < 0.7. The complete selection has an efficiency of about 
58 % for hadronic Z° decays, and selected a total of 1 923 240 events in the data. Of these events, 

1 A right handed coordinate system is used, with positive z along the e~ beam direction and x pointing towards the 
centre of the LEP ring. The polar and azimuthal angles are denoted by 8 and (j>, and the origin is taken to be the centre 
of the detector. 
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4.9% were recorded below the peak of the Z° resonance at a centre-of-mass energy of E cm = 89.4 GeV 
and 7.4% above the peak at E cm = 93.0 GeV. Calculations using ZFITTER || indicate that these 
off-peak events change the measured value of i?b by —0.00004, a correction which is —0.02 % of the 
measured value and is not applied to the result presented here. 

The event selection is designed to have the same efficiency for each quark flavour. However, bb 
events have a higher average charged particle multiplicity than other flavours, and hence a slightly 
higher probability of passing the event selection requirement of at least seven charged tracks. Owing 
to the high mass of the b quark, bb events also have a slightly different thrust distribution, but Monte 
Carlo studies showed this to have a much smaller effect than that from the track multiplicity cut. 
According to the Monte Carlo, these biases increase the measured value of i?b i n the selected event 
sample by 0.32 % of its value. Non-hadronic background (mainly e + e~ — * t + t~ events) constitutes 
about 0.065 % of the selected events, and reduces the measured value of by 0.065 %. The combi- 
nation of these event selection effects increases the measured value of R\, by (0.25 ± 0.15)%, where 
the error is due to the uncertainty in the simulation of the track multiplicity distributions and Monte 
Carlo statistics 

Charged tracks and electromagnetic calorimeter clusters with no associated track were combined 
into jets using a cone algorithm with a cone half angle of R = 0.55 rad and a minimum jet energy 
of 5 GeV. This algorithm, rather than the JADE algorithm employed in Q, increases the fraction of 
tracks in the jet coming from the b hadron decay, and improves the b-tagging performance. 

Monte Carlo simulated events were used for evaluating backgrounds, acceptances for charm and 
light quark events, and efficiency correlations between the two hemispheres of an event. Hadronic 
events were simulated with the JETSET 7.4 generator ||, with parameters tuned by OPAL Q. The 
fragmentation function of Peterson et al. was used to describe the fragmentation of b and c 
quarks. The generated events were passed through a program that simulated the response of the 
OPAL detector (T^j and the same reconstruction algorithms as the data. Further detailed studies of 
the tracking performance in the data were made, and the results were used to tune the simulated 
tracking performance. 

4 Vertex Tagging 

Hadronic Z° decays into b quarks are tagged by reconstructing secondary vertices significantly sepa- 
rated from the primary vertex, taking advantage of the relatively long (~ 1.5 ps) lifetime, high decay 
charged multiplicity and high mass (~ 5 GeV) of b hadronsf]. Information characterising these fea- 
tures of b hadron production and decay is combined using an artificial neural network algorithm to 
produce a single vertex tagging variable for each hemisphere. 

4.1 Primary Vertex Reconstruction 

Primary vertices were reconstructed for each event using the algorithm described in but applied 
separately to each hemisphere of the event. Reconstructing a separate independent primary vertex in 
each hemisphere avoids an efficiency correlation from sharing a single primary vertex reconstructed 
using tracks from both hemispheres |l^]. The position and uncertainty of the beam spot (e + e~ collision 
region) was used as a common constraint in both hemispheres, but its effect on the correlation is 
small. The beam spot position was measured using charged tracks from many consecutive events, 
thus following any significant shifts in beam position during a LEP fill 

All good tracks in one hemisphere were fitted to a common vertex in three dimensions. The tracks 
were required to satisfy various requirements to ensure they were well measured, including having at 
least 20 hits in the jet chamber, a transverse momentum with respect to the beam axis of at least 
0.15 GeV and a distance of closest approach to the beam spot in the r-<p plane of less than 5 cm. Tracks 

2 The convention c = 1 is assumed throughout this paper. 
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with a large \ 2 contribution to the primary vertex fit were removed one by one until each remaining 
track contributed less than 4 to the \ 2 ■ I n about 0.7% of the events, no tracks remained after this 
procedure, in which case the beam spot position was used. 

Using the hemisphere primary vertex reconstruction, the primary vertex resolution in Monte Carlo 
simulated events with the 1994 silicon detector is about 70 ^m in x, 15 /iin in y (where it is dominated 
by the beam spot constraint) and 110 ^m in z. Using the event primary vertex reconstruction, the 
resolutions in x, y and z are 50 /im, 15 fim and 80 fim respectively. The poorer primary vertex resolution 
in the hemisphere reconstruction method leads to a non-b impurity that is about 10 % higher (relative) 
for a given b-tagging efficiency, but a much reduced hemisphere efficiency correlation and associated 
systematic error. 



4.2 Secondary Vertex Reconstruction 

A three dimensional version of the algorithm described in jij was used to reconstruct a secondary 
vertex in each jet. All tracks in the jet with momentum p > 0.5 GeV, a distance of closest approach 
to the hemisphere primary vertex in the r-cft plane ('impact parameter') do < 0.3 cm, and an impact 
parameter error a^o < 0.1 cm, were fitted to a common vertex in three dimensions. These requirements 
preferentially select vertices from b hadron decays because of the high decay multiplicity and hard 
fragmentation of the b quark. Tracks with a large x 2 contribution to the fit were removed one by 
one until each track contributed less than 4 to the overall \ 2 ■ At least three tracks were required to 
remain in the fit for the secondary vertex finding to be considered successful. 

The vertex decay length L was calculated as the length of the vector from the primary to the 
secondary vertex, the vector being constrained to lie along the direction of the jet axis. L was given 
a positive sign if the secondary vertex was displaced from the primary in the direction of the jet 
momentum, and a negative sign otherwise. Vertices with L > ('forward') were used in the tagging 
of b hadron decays, and those with L < ('backward') were used to reduce the systematics associated 
with the detector resolution. 



4.3 Vertex Tagging Neural Network 

Jets containing a vertex with at least three tracks and satisfying \L/ai\ > 3 (where ctl is the estimated 
error on L, derived from the track error matrices) were considered pre-selected and input to a neural 
network algorithm [14]. In the Monte Carlo, the efficiency of this pre-selection for b-jets is about 49% 
and the b-purity is about 76%. 

The neural network was used to enhance the b-purity, and has five input nodes, eight hidden nodes 
and one output node. The first three inputs were derived directly from the reconstructed secondary 
vertex: the decay length significance L/ai, the decay length L, and the number of tracks in the 
secondary vertex N s . The fourth input to the network tests the stability of the vertex against mis- 
measured tracks. The secondary vertex track with the highest impact parameter significance with 
respect to the primary vertex (do/cr^) was removed from the secondary vertex, and the secondary 
vertex fit repeated resulting in the 'reduced' decay length significance Lr/(Jl r - For genuine b hadron 
decays, Lr/(Jl, r is large, whilst for vertices caused by one high impact parameter mis-measured track 
(which is removed from the 'reduced' vertex), Lr/gl r is small. 

The fifth input exploits the high mass of b hadrons compared with charm hadrons, using a method 



similar to that described in [12]. For each track in the jet, a weight X that it came from a b hadron 
decay was computed, using a separate artificial neural network trained on bb Monte Carlo events. This 
network has six inputs, eight hidden nodes and one output node. The inputs are: the scaled track 
momentum x p = p/i^eam! the sine of the angle of the track with respect to the jet axis sin 6^ = pt/p; 
the impact parameter significances of the track with respect to the reconstructed secondary vertex in 
the r-(p ((do /& d ) sec) and r-z ((zo/a Zo ) sec ) planes; and the impact parameter significances of the track 
with respect to the hemisphere primary vertex in the r-</> ((do/cd )prim) and r-z ((-2o/cz )prim) planes. 
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The network was trained separately on tracks with no silicon hits, silicon r-<j> hits only, and both silicon 
r-(j) and r-z hits. The output X of this network peaks close to zero for tracks from fragmentation, and 
close to one for tracks from b hadron decay. 

All tracks within a jet were then ordered by decreasing X [i.e. most b hadron decay-like tracks 
first), and the first two were clustered together. Other tracks were added in turn to the cluster, until 
the invariant mass of the tracks (assuming them to be pions) exceeded the charm hadron mass, taken 
to be m_o = 1.9 GeV. The value Xrj of X for the track which caused the cluster invariant mass to 
exceed was then used as the fifth input for the vertex neural network. For charm hadron decays, 
Xd is usually small, since tracks from fragmentation have to be included to exceed mo- In contrast, b 
hadrons have enough mass that the threshold can usually be exceeded with tracks from the b hadron 
decay alone, leading to a value of Xr> close to one. If the cluster mass did not exceed after all 
tracks in the jet had been added, Xp was set to zero; this happens more often in charm and light 
quark jets than in b jets. 

4.4 Vertex Tag Definition 

Samples of five-flavour Monte Carlo jets passing the pre-selection and with L > were used to train 
the main neural network (with inputs L/ctl, L, N s , Lr/ctl r and Xrj), to produce output distributions 
peaking close to zero for light flavour and charm jets, and outputs close to one for b jets. 

In order to reduce sensitivity to the modelling of the detector resolution, the technique of 'folding', 
i.e. subtracting the number of hemispheres tagged with negative L, was employed H. This technique 
works well if the tagging variable (in this case the vertex tag neural network output) is symmetric 
about zero for jets containing no particles with detectable lifetime. To achieve this symmetry, the 
variables L/ctl, L and Lr/(jl r were modified before being input to the neural network. For L/gl and 
L , the absolute values \L/ai\ and \L\ were taken, while Lr/<jl r was signed positive if it originally 
had the same sign as L, and negative otherwise. The impact parameter significances used to calculate 
Xu also had their signs reversed if L was negative. The magnitude of the vertex tagging variable B 
was then defined as \B\ = — ln(l — 6) where 6 is the raw neural network output (between zero and 
one), and the sign of B was taken to be the sign of L. The logarithmic transformation is used to 
expand the scale of the tagging variable in the region just below 6 = 1. 

The use of folding also requires that equations [I] and [2] are modified appropriately. The tagging 
efficiencies e are replaced by the difference of forward and backward tagging efficiencies e v — e-. The 
number of tagged hemispheres Nt is replaced by the difference between the numbers of forward and 
backward tagged hemispheres N v — N-, and the number of double tagged events N tt is replaced by 
N vv — N V y + N-^- where N vv , N Y ^ and Nyy are the numbers of events with two forward tags, one 
forward and one backward tag, and two backward tags respectively [Q]. 

As the silicon detector did not provide r-z information in 1992, a separate version of the vertex 
tagging algorithm was used for these data. Primary and secondary vertices were reconstructed in 
the r-(f) plane only, and the r-z impact parameter significances were not used in the calculation of 
Xd. According to the Monte Carlo, the tagging impurity for 1992 data using the r-<j) only tagging 
algorithm is about 30 % higher than that for 1994 data at the same b-tagging efficiency. 

4.5 Vertex Tag Performance 

The distributions of the input variables for the track neural network used to derive Xp are shown in 
Figure |], and the distributions of the vertex tag neural network input variables are shown in Figure |^, 
for 1994 data and Monte Carlo simulation. In general the modelling of the input variables is good, but 
some discrepancies are visible, particularly in the distributions of N s and Xjj. These are attributed 
to imperfections in the detector modelling and in the simulation of b hadron decays. The effects of 
the former are addressed in the systematic errors, whilst the latter have no impact on the result since 
the b-tagging efficiency is determined from the data and not from the Monte Carlo. 
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Fi gure 1: Distributions of the track neural network input variables Xp, sin^, (<io/<7<io)sec> {zo/cr zo ^ sec , 
(^o/cdo)prim and (zo/a ZQ ) pT i m used to derive Xp- The sum of the distributions for the three track 
classes (no silicon hits, silicon r-(f) hits only, and both silicon r-<p and r-z hits) are shown for 1994 data 
(points) and 1994 Monte Carlo simulation (histograms). The estimated contribution of tracks from b 
decay, b fragmentation, and charm and light quark background are indicated. 
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Figure 2: Distributions of the vertex tagging neural network inputs L/ctl, L, Lji/aL R , N s and Xjj for 
1994 data (points) and 1994 Monte Carlo simulation (histogram). The contributions from uds, c and 
b jets are indicated. 
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The distribution of the vertex tagging variable B is shown for 1994 data and Monte Carlo in 
Figure ||. Good qualitative agreement is seen. The Monte Carlo distribution for uds jets is seen to 
be approximately symmetric about B = 0, as required for the folding procedure to work well. The 
asymmetry in the uds distribution is caused by jets containing a long-lived strange particle (Kg, A or 
other hyperon) and by uds jets containing a gluon splitting to a bb or cc pair. The distributions for 
charm and especially for b jets are shifted to large positive values of B, and at high values the b jet 
purity exceeds 99%. 

The hemisphere vertex tag was defined from the vertex tag of any jet in the hemisphere. If more 
than one jet in the hemisphere passed the vertex tag preselection, the one with the highest value of 
\B\ was used. A hemisphere was then defined to be forward tagged with the vertex tag V if B > 2.8 
and backward tagged (V) if B < —2.8. The value of 2.8 was chosen to minimise the overall error on 
i?b when also including the lepton tag. 

5 Lepton Tagging 

Leptons with high momentum p, and a large momentum component transverse to the jet axis pt, are 
expected to come mainly from semileptonic decays of b hadrons, because of the hard fragmentation 
and high mass of the b quark. Electron candidates with p > 2 GeV and pt > 1.1 GeV, and muon 
candidates with p > 3 GeV and pt > 1.4 GeV, were used to tag bb events. Both electrons and muons 
were restricted to the polar angle range | cos 6*| < 0.8 to ensure a well-defined acceptance and good 
Monte Carlo modelling of the efficiencies and backgrounds. 

The sources of lepton candidates are divided into two classes. The first class consists of prompt 
leptons from the decays of b and c hadrons (including b — > r — > i and b — > J/ip — > t and those from b 
and c hadrons produced in gluon splitting). These leptons are included in the definitions of e b , e c and 
£ uds £ Qr i e pt on tag. Monte Carlo is used to estimate e° and e uds , whereas e b is determined from the 
data. The second class of lepton candidates consists of everything else: real leptons produced from 
Dalitz decays, photon conversions and the decays in flight of K 1 * 1 and ir^, and hadrons mis-identified as 
leptons. This background is estimated using a combination of data and Monte Carlo, and is subtracted 
from the number of identified lepton candidates before input to the fit for R^. These lepton candidates 
are therefore not included in the definitions of e b , e c and e uds . 

5.1 Electron Identification 

Electrons were identified using an artificial neural network. The algorithm is a simplified version of 
that described in 113], using only six rather than twelve neural network inputs. The inputs are: the 
momentum and polar angle of the track, the energy-momentum ratio E/p, the number of electromag- 
netic calorimeter blocks contributing to the energy measurement, the normalised ionisation energy loss 
^ae/Ax an d the error on the ionisation energy loss cr dE / dx . The normalised dE/dx value is defined 
as N^ E j Ax = (dE/dx — (dE /dx)o) / (a dE / dx )o, where (dE/dx)o is the value and (o"d£;/dx)o the error 
expected for the track, assuming it to be an electron of the measured track momentum. The neural 
network output for electron candidates in data and Monte Carlo is shown in Figure 0(a). Electron 
candidates were required to have a neural network output greater than 0.95. 

Photon conversion candidates were rejected using another neural network algorithm, using spatial 
matching, invariant mass and momentum information of the electron candidate and an oppositely 
charged partner track. This algorithm is similar to that described in |15| |, but uses the new electron 
neural network algorithm described above. The distributions of this neural network output are shown 
in Figure ||(b). Electron candidates were required to have a conversion tagging neural network output 
smaller than 0.5. After all these requirements, the expected identification efficiency for electrons from 
decays of b hadrons within the kinematic and geometrical acceptance is about 68 %. 

The non-prompt background in the tagged electron sample consists of hadrons mis-identified as 
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tagging variable B 



Figure 3: Distributions of the vertex tagging variable B for the 1994 data (points) and 1994 Monte 
Carlo simulation (histogram). The positions of the forward and backward tag cuts are shown by 
dashed lines. The contributions from uds, c and b jets are indicated. 
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Figure 4: Performance of neural network electron identification in data and Monte Carlo: (a) Nor- 
malised distributions of electron neural network output iV c i for all tracks with N^ E ^ dx > —2, p > 2 GeV 
and pt > 1.2 GeV; (b) Normalised distributions of the photon conversion tagging neural network out- 
put for identified electron tracks with N e \ > 0.95. In each case, the 1994 data is shown by the points 
with error bars, and the Monte Carlo contributions from prompt electrons, photon conversions and 
mis-identified hadrons by the histograms. The selected regions are shown by the dashed lines and 
arrows. 
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1123 


1329 
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1289 






-AW 


+ AW 


4162 


4993 


10101 


5369 



Table 1: Numbers of hadronic events, tagged hemispheres and double tagged events selected in each 
year of the data. Background has been subtracted from the lepton samples, and the resulting counts 
are quoted to the nearest whole number. 



electrons, untagged photon conversions and a small number of electrons from Dalitz decays of light 
mesons. The background mis-identification probabilities were found to depend strongly on the p and 
Pt of the track, but were otherwise largely independent of the event flavour. Since the track p and 
Pt distributions are different in bottom, charm and light quark events, the number of fake leptons in 
hemispheres opposite tagged and untagged hemispheres are different. Therefore, the Monte Carlo was 
used to determine the fake probability per track as a function of p and pt- The probabilities were 
combined with the two dimensional distributions of track p and pt measured opposite untagged and 
tagged hemispheres in the data, to estimate the number of fake electrons expected opposite each type 
of hemisphere. These estimates were then subtracted from the number of lepton tagged hemispheres 
and events before input to the fit for R^. 

In total, 26185 hemispheres were tagged by electrons after the photon conversion rejection. Of 
these, 918±30 were attributed to hadronic fakes, 784±28 to untagged photon conversions and 210±14 
to Dalitz decays of light mesons, where the background rates have been evaluated by subdividing the 
fake probabilities by source as a function of p and pt- The errors are due to data statistics only. 



5.2 Muon Identification 

Muon candidates were identified by matching track segments reconstructed in the four-layer external 
muon chambers to tracks extrapolated from the central tracking detectors. The measured dE/dx was 



also required to be consistent with a muon. The algorithm is described in detail in [16]. The expected 
identification efficiency for muons from decays of b hadrons within the kinematic and geometrical 
acceptance is about 79 %. 

The muon background was estimated from the Monte Carlo using the same techniques as described 



in Section 5.1. In total, 21558 hemispheres containing muon candidates were found in the data, of 



which 3311 ± 58 were attributed to hadronic background. The errors are due to data statistics only. 



6 Measurement of 

The numbers of hadronic events, tagged hemispheres and double tagged events found in each year 
of the data are listed in Table |lj. The symbol iVj represents the number of hemispheres tagged by 
tag i, where i = v for forward tagged vertices, i = v for backward vertices, i = i for leptons and 
i = a for hemispheres tagged by either a forward vertex or a lepton. The symbol represents the 
number of events tagged by tag i in one hemisphere and tag j in the other hemisphere. The expected 
photon conversion and hadronic backgrounds have already been subtracted from the lepton counts as 
described in Section ||. 

The hemisphere tagging probabilities for charm and light quark events, estimated using Monte 
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Tag 


Year 


cc(%) 


uu±dd±ss(%) 


Vertex 


1992 


0.463 ± 0.009 


0.0279 ± 0.0034 




1993 


0.457 ±0.015 


0.0256 ± 0.0020 




1994 


0.473 ± 0.007 


0.0256 ± 0.0014 




1995 


0.491 ± 0.009 


0.0253 ± 0.0023 


Lepton 


1992 


0.425 ± 0.009 


0.0181 ±0.0016 




1993 


0.455 ± 0.013 


0.0196 ±0.0010 




1994 


0.443 ± 0.006 


0.0172 ± 0.0006 




1995 


0.445 ± 0.008 


0.0155 ±0.0010 


Combined 


1992 


0.890 ±0.013 


0.0467 ± 0.0037 




1993 


0.902 ±0.016 


0.0425 ± 0.0030 




1994 


0.918 ± 0.009 


0.0425 ± 0.0015 




1995 


0.936 ±0.012 


0.0405 ± 0.0025 



Table 2: Hemisphere tagging probabilities for charm and light quark events estimated from the Monte 
Carlo for each year of the data. The errors are due to Monte Carlo statistics only. 



C b - 1 (%) 


Vertex 


Lepton 


Combined 


1992 


0.84 ±0.24 


1.82 ± 1.29 


1.01 ±0.21 


1993 


1.10 ±0.26 


1.34 ± 1.28 


1.13 ±0.20 


1994 


0.77 ±0.22 


1.13 ± 1.25 


0.92 ±0.18 


1995 


0.99 ±0.24 


0.92 ± 1.28 


0.95 ±0.19 



Table 3: Hemisphere efficiency correlations in bb events estimated from the Monte Carlo with small 
corrections for each year of the data. The errors are due to Monte Carlo and data statistics only. 



Carlo simulation, are given in Table ^. The tag probabilities vary slightly from year to year because of 
the differences in silicon geometrical acceptance and the r-(f> only tag used in 1992. The Monte Carlo 
predicted no significant difference in tagging efficiency between uu, dd and ss events, and the effect 
on i?b of assuming a common efficiency was estimated to be less than 10 -6 . 

The bb tagging efficiency correlations C b — 1 for each year of the data are listed in Table |3|. These 
values were determined from a Monte Carlo sample of 7 million bb events, by calculating the ratio 
of the event double tagging probability to the product of the hemisphere single tagging probabilities. 
Small corrections for the differing detector performance in each year were derived from the data, as 



will be discussed in Section p.5| . The errors include contributions from the Monte Carlo and data 
statistics only. 

The values of and the tagging efficiencies e b derived for each year of the data are given in 
Table The results are obtained by solving equations |l| and The full result from the combined 
tag is calculated using the number of tagged hemispheres N a — Ny and the number of double tagged 
events iV aa — iV av - ± Nyy. Results are also given for the vertex tag alone (N v — Ny and iV vv — N v - + N—) 
and the lepton tag alone (Ng and iV«). The results have been corrected for the event selection bias 
described in Section ||. Only the data statistical errors are included. The results from the combined 
tag for the four years agree with each other at a x 2 °f 1-8 for 3 degrees of freedom; the corresponding 
values of \ 2 f° r the vertex and lepton tags alone are 3.4 and 2.9 respectively. Combining the data 
from all four years, the value of is measured to be: 

R h = 0.2178 ±0.0011 

where the error is due to the data statistics only. 
The result depends on R c as follows: 



0.056- 



Rb R, 



c 
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Vertex 


Lept 


on 


Combined 


.Rb 1992 


n 91 fio 


it 


U.UUOO 


n 91 ^fi 


_i_ 
rt 


0.0174 


n 91 fi/i 


j_ 


fl DD97 
U.UUZ / 


1993 


fl 91 89 


_i_ 
in 


n nn^n 

U.UUoU 




_i_ 


0.0199 


fl 91 89 


_1_ 


U.UUZD 


1994 


0.2176 


± 


0.0021 


0.2113 


± 


0.0110 


0.2194 


± 


0.0018 


1995 


0.2120 


± 


0.0027 


0.2171 


± 


0.0160 


0.2157 


± 


0.0023 


Combined 


0.2163 


± 


0.0013 


0.2184 


± 


0.0074 


0.2178 


± 


0.0011 


e b 1992 


0.1893 


± 


0.0030 


0.0439 


± 


0.0035 


0.2254 


± 


0.0029 


1993 


0.2002 


± 


0.0029 


0.0406 


± 


0.0032 


0.2369 


± 


0.0028 


1994 


0.2087 


± 


0.0021 


0.0477 


± 


0.0025 


0.2429 


± 


0.0020 


1995 


0.2211 


± 


0.0030 


0.0477 


± 


0.0035 


0.2550 


± 


0.0028 



Table 4: Values of Rb and e b measured in each year of the data, after correlation and event selection 
bias correction. Only statistical errors are included. 



where AR C is the deviation of R c from the value 0.172 predicted by the Standard Model and used in 
this analysis. 

The systematic errors coming from sources other than R c are discussed below and are summarised 
in Table ||. Most of the systematic errors arise through the charm and light quark hemisphere tagging 
efficiencies, e c and e uds , and through the bb tagging efficiency correlation C b . The dependence of i?b 
on these quantities is given approximately by: 

Ai? b Ae c Ae uds AC b 

— * = _ .059— - 0.010 — g- + — j- (3 

R h e c e uds C h 

The systematic errors arising from the charm and light quark efficiencies are discussed in detail in Sec- 
tion 0, and those from the efficiency correlation in Section ||. The systematic errors on the efficiencies 
are also given in Table |5|. The only other source of systematic error is the hadronic event selection, 
which was discussed in Section ||] and gives rise to an error of ±0.00033 on R^. 

As a cross check, the cut on the vertex tag neural network output B was varied in the range 2.2-3.8, 
and the cuts on the lepton transverse momenta were varied by up to ±0.3 GeV from their nominal 
values. The values of R\> obtained are shown in Figure ^ together with the uncorrelated parts of the 
statistical and systematic errors. No significant trend in the measured value of i?b is observed within 
these errors. 



7 Systematic errors: tagging efficiencies 

Systematic errors on the tagging efficiencies e c and e uds arise from the understanding of the tracking 
detectors, the electron and muon identification, and the various physics parameters input to the Monte 
Carlo simulation. 



7.1 Tracking detector performance 

The evaluation of the charm and light quark tagging efficiencies for the vertex tag requires an accurate 
simulation of the detector resolution for charged tracks. The Monte Carlo simulation has been tuned 
to reproduce the tracking resolution seen in each year of the data by studying the impact parameter 
distributions of tracks, as functions of track momentum, polar angle and the different sub-detectors 
contributing hits. 

The effect of uncertainties in this procedure was evaluated as follows: 

Tracking resolution: The sensitivity to the tracking resolution was assessed by degrading or improv- 
ing the resolution in the Monte Carlo. This was done by applying a single multiplicative scale 
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Source 


Ae c /e c (%) 


A uds / uds ty\ 


Ai? b 


Tracking resolution 


1.24 


4.0 


0.00017 


Tracking efficiency 


0.80 


4.0 


0.00014 


Silicon hit matching efficiency 


0.82 


2.8 


0.00009 


Silicon alignment 


0.58 


2.1 


0.00008 


Electron identification efficiency 


1.11 


0.5 


0.00015 


Muon identification efficiency 


0.64 


0.2 


0.00009 


c quark fragmentation 


2.26 




0.00028 


c hadron production fractions 


3.66 




0.00046 


c hadron lifetimes 


0.55 




0.00007 


c charged decay multiplicity 


1.09 




0.00014 


c neutral decay multiplicity 


2.39 




0.00030 


Branching fraction B(D — > K°) 


1.20 




0.00015 


c semileptonic branching fraction 


2.44 




0.00031 


c semileptonic decay modelling 


2.34 




0.00029 


Gluon splitting to cc 


0.34 


6.3 


0.00018 


Gluon splitting to bb 


0.50 


9.3 


0.00027 


K° and hyperon production 




0.3 


0.00001 


Monte Carlo statistics (c, uds) 


0.66 


2.5 


0.00010 


Subtotal Ae c and Ae uds 


6.65 


13.3 


0.00090 


Electron identification backgrounc 






0.00039 


Muon identification background 






0.00041 


Efficiency correlation AC b 






0.00066 


Event selection bias 






0.00033 


Total 


0.00129 



Table 5: Systematic errors on the measured value of R^. The uncertainties on the charm and 
light quark efficiencies for each tag are also given. The systematic errors arising from the efficiencies 
and lepton identification background are discussed in Section [7], those from efficiency correlation in 
Section || and that from the event selection in Section |3|. 



factor /3 to the difference between the reconstructed and true track parameters of all charged 
tracks jl. A ±10% variation was applied to the r-(j) track parameters (the impact parameter do 
and track azimuthal angle (fio). An independent ±10% variation was applied to the analogous 
parameters in the r-z plane. Together with uncertainties in the simulation of b hadron decays 
and fragmentation, these variations can account for all the discrepancies observed between data 
and Monte Carlo in the neural network input and output distributions. The result of these 
variations is an error of ±0.00017 on R^. 

Track reconstruction efficiency: The reconstruction efficiency for charged tracks is estimated to 



exceed 98% for the momentum and impact parameter requirements described in sections 4.1 



and |4.2 . Most of the losses occur in small regions of <fi around the jet chamber cathode and 
anode wire planes. The <f> distribution of reconstructed tracks is well reproduced by the Monte 
Carlo simulation, as is the two-track resolution and the rate of tracks lost because they lie on 
top of another track when reflected in the anode plane. Residual discrepancies indicate that the 
track reconstruction efficiency is modelled to within ±1 %. The effect on R^ was assessed by 
randomly discarding 1 % of tracks in the Monte Carlo, and was found to be ±0.00014. 

Silicon hit association efficiency: The Monte Carlo was tuned to model the overall rate of associ- 
ating silicon hits to tracks, and to model known inefficient regions of the silicon detector in the 
1993 and 1994 data. The residual discrepancies between data and Monte Carlo association effi- 
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Figure 5: Values of R^, obtained at different cut values. The central value is indicated by the arrow 
and the dashed lines indicate the total statistical and systematic error. The error bars indicate 
the uncertainties on the differences from the central result, including both statistical and systematic 
components. 



ciencies were found to be within ±1 % for the r-<p and ±3% for the r-z hits. The hit association 
efficiency in the Monte Carlo was varied within these errors resulting in a systematic error on 
R h of ±0.00009. 

Silicon alignment: The position of the r-(f> silicon wafers in the azimuthal direction, and the r-z 
wafers along the z direction, is determined by an alignment procedure using Z° — > events 
to a precision of about 10 /um ||. This uncertainty is included in the Monte Carlo simulation, 
but has little impact on the track resolution in hadronic events which is dominated by multiple 
scattering. The radial alignment uncertainty is much more important since radial shifts of 
individual wafers can lead to systematic mis-measurement of the decay length L in jets contained 
largely within a single wafer. The radial alignment was studied using cosmic rays (which are 
incident on the wafers at all angles) recorded throughout the data taking period and found 
to be good to a precision of ±20 fim. The effect on the tagging efficiencies was studied by 
systematically displacing one or both silicon barrels radially by 20 /im in the simulation, and 
was found to correspond to an error of ±0.00008 on i?b- 

The forward, backward and folded vertex tagging rates measured in each year of the data are shown 
in Figure ||. The Monte Carlo prediction as a function of the r-<j) resolution parameter [5 is also shown. 
The Monte Carlo forward tagging rate is not important, since the b-tagging efficiency is determined 
directly from the data using the double tagging technique. It is always somewhat higher in the Monte 
Carlo than in the data, but the observed discrepancies are within the uncertainties associated with 
b-physics and detector simulation in the Monte Carlo. The changes in tagging rates as a function 
of year are mainly caused by differences in the silicon detector in each year of data taking. In 1992, 
only r-cj) coordinate information was available, the r-z wafers being installed before the 1993 run |J. 
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Figure 6: Forward (/ v ), backward (f-) and folded (/ v — f-) hemisphere tagging fractions in each 
year for the vertex tag. The data tagging rates are shown by the shaded bars, the width of the bars 
representing the statistical error. The Monte Carlo predictions as a function of the r-<p resolution 
scaling parameter /3 r9 i are shown by the points with error bars. 



18 



A faulty silicon module in 1993 was replaced for the 1994 run, and a new detector with increased 
coverage installed for 1995. 

The Monte Carlo simulation has been tuned by studying the impact parameter distributions of 
single tracks. This procedure is not sensitive to possible coherent effects affecting all the tracks in 
a small number of jets in a correlated way. For example, such effects can be caused by a badly 
mis-reconstructed primary vertex, which could change the impact parameters of all tracks in a jet, 
potentially giving an apparent large negative decay length and hence a backward tag. However, the 
discrepancies in the backward tagging rates seen in Figure || are within the range covered by the 
±10% variation in the r-(j) resolution parameter f3, so such coherent tracking resolution effects are not 
likely to be important. This was also checked by studying tails in the distributions of differences in 
the primary vertex position measured in each hemisphere. 

The effect of the r-</) resolution scaling on the light quark and charm tagging efficiencies is shown 
in Figure]^, for variations in f3 r ^ between 0.8 and 1.2 (±20%). These variations cause large changes 
in the light quark (uds) forward and backward tagging efficiencies, but these changes largely cancel 
in the folded tagging efficiencies, thus reducing the systematic uncertainty. In contrast, the charm 
tagging efficiency is dominated by real lifetime tags, and the folding procedure has only a small effect. 



7.2 Electron identification 

Monte Carlo simulation was used to predict the efficiency for identifying prompt electrons from c — > e 
decays and those from gluon splitting g — ► (bb, cc) — > e, and background rates as function of track 
p and pt- The modelling of the electron identification requirements in the Monte Carlo was checked 
by studying in detail the distributions of the input variables of the neural network. The effects on 
the electron identification performance of discrepancies due to each input variable were added in 
quadrature to estimate the total systematic errors. This procedure was cross checked with samples of 
electrons and hadrons selected from various pure control samples. 

The Monte Carlo efficiency prediction was checked by studying the neural network input distribu- 
tions and with samples of pure electrons from e + e~ — ► e + e~ events and photon conversions, and was 
found to be modelled to a relative precision of ±4%. According to the Monte Carlo, about 2.6% of 
true prompt electrons are rejected by the photon conversion finder, and this was checked to a relative 
precision of ±30% by comparing samples of high p and pt electrons in the data and Monte Carlo. 
This corresponds to a further error of ±0.8 % on the prompt electron efficiency. Combining both the 
electron identification and photon rejection uncertainties gives a systematic error on of ±0.00015. 

The two main non-prompt backgrounds in the tagged electron sample are charged hadrons (mainly 
ir^ 1 ) which are mis-identified as electrons, and untagged photon conversions. The effect of mis- 
modelling of the neural network input distributions in the Monte Carlo was found to correspond 
to an uncertainty of ±21 % in the hadron mis-identification probability. Consistent results were found 
using control samples of selected in Kg —* 7t + tt~ and r — > 37r decays. The probability for mis- 
identified hadrons to pass the conversion rejection requirements was found to be modelled to ±5 % 
using the Kg — ► 7r + 7r~ sample, a sample of prompt muons (which have similar p and pt distributions 
to prompt electrons and no photon conversion background) and an inclusive electron-depleted sample 
selected using an anti-cut on the electron identification neural network. 

The untagged photon conversion background was studied using a sample of identified electrons, 
enriched in conversions by requiring the tracks to have low p and pt- The number of untagged 
conversions was found to be modelled to a precision of ±15%. A small number of electrons are also 
produced in the Dalitz decay of 7r° and rj mesons. The Monte Carlo prediction was used for this rate, 
and a ±20 % uncertainty used to assess the systematic error, based on measurements of inclusive it 



and rj production in hadronic events [17|. The total systematic error on from all the background 
sources was found to be ±0.00039. 
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Figure 7: Forward (e), backward (e) and folded (e — e) vertex tagging efficiencies, and their differences 
for light quark (uds) and charm hemispheres in f 994 Monte Carlo, as a function of the r-(p resolution 
scaling parameter (3^. The shaded bands around the central f3 r ^, = 1.0 points show the total statistical 
error, and the error bars on the other points show the statistical errors on the difference between them 
and the central points. 
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7.3 Muon identification 



The muon identification efficiency was studied using various control samples in data and Monte Carlo 
[|l~6|]. The muon matching requirements were studied using muon pairs from two-photon production, 
for muons in the range 2 to 6 GeV, and muon pairs from Z° —* fJ- + fJ.~ events for muons above 30 GeV. 
The efficiency of the dE/dx requirement was studied using Kj? — » it + it~ decays. The Monte Carlo was 
found to model the matching and dE/dx requirement efficiencies to relative precisions of 2.1% and 
2.2 % respectively, to give a total error on the muon identification efficiency of 3%. This corresponds 
to an error of ±0.00009 on R^. 

The various sources of background muons were studied using control samples from Kg —* 7r + 7r~ 



and r — > 3ir decays, as in [16]. The background in the data was found to be a factor 1.13 ±0.09 higher 
than that in the Monte Carlo, and the fake probabilities per track were corrected accordingly. The 
resulting uncertainty on i?b is ±0.00041. 

7.4 Simulation input parameters 

The charm and light quark efficiencies are sensitive to the following simulation input parameters: 

Charm quark fragmentation: The charm tagging efficiency e c increases with the scaled energy 
of the weakly decaying charm hadron. The mean scaled energy (xe) of charm hadrons 



produced in cc events at LEP has been measured by ALEPH, DELPHI and OPAL [18 



These measurements have been averaged by the LEP electroweak working group to give a value 
of (xe) = 0.484 ± 0.008 for weakly decaying charm hadrons [21]. The effect of this uncertainty 
was assessed by reweighting events in the Monte Carlo changing (xe) within this range. The 
fragmentation functions of Peterson J1Q|] , Collins and Spiller ]|22| , Kartvelishvili |23| and the 
Lund group [24 were used as models to determine these event weights. The largest variation 
was found using the model of Collins and Spiller, and the resulting variation in of ±0.00028 
was assigned as a systematic error. The charm fragmentation also affects the tagging efficiency 
via the number of tracks produced in the fragmentation process, but this effect was found to be 
much smaller than the direct energy dependence and was neglected. 

Charm hadron production fractions: Because of the different charm hadron lifetimes and decay 
modes, the vertex tagging efficiency in cc events depends on the mixture of weakly decaying 
charm hadrons. The tagging efficiency for D + mesons is approximately three times that for D°, 
whilst that for D+ mesons is approximately 15% higher than D° and that for A+ is only 15% 
of that for D°. The fractions of D + , D°, D+ and A+ were varied according to the production 
fractions measured at LEP [19,25| as averaged by the LEP electroweak working group [21]. The 



contribution from A+ was scaled by 1.15 ± 0.05 to account for other weakly decaying charm 
baryons. The errors were combined taking their correlations into account to give an error on 
i?b of ±0.00031. The dependence of the tagging efficiency on the fraction of weakly decaying 
charm hadrons produced via the decay of excited charm states (D* and D**) was found to be 
negligible. 

The systematic error on the lepton tagging efficiency in cc events is derived from the inclusive 
charm semileptonic branching fraction, which is in turn derived from the individual charm hadron 
semileptonic branching fractions (see below). This introduces an additional dependence on the 
charm hadron production fractions which is correlated with that from the vertex tag, taking the 
total error on i?b from this source to ±0.00046. 

Charm hadron lifetimes: The lifetimes of the weakly decaying charm hadrons were varied sepa- 
rately within the errors quoted by the Particle Data Group [26|. Their contributions to the error 



on i?b were added in quadrature to give an error on i?b of ±0.00007. 
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Charm hadron charged decay multiplicity: The reconstructed secondary vertex track multiplic- 
ity N s is required to be at least three, and is also used as an input to the vertex tag neural 
network. The tagging efficiency in cc events therefore depends strongly on the charged track 
multiplicity of charm hadron decays. The average charged track multiplicity of D + , D° and D+ 
decays (including the charged decay products of any produced K° mesons) has been measured 
by MARK III p7| . The average multiplicity in Monte Carlo events was varied within the range 
given by MARK III using several different reweighting schemes, keeping the inclusive branching 
ratios to K° and A constant for each charm hadron. For charm baryons, for which no measure- 
ments are available, a variation of ±0.5 was taken. The variations for each charm hadron species 
were combined in quadrature to give an overall error on of ±0.00014. 

Charm hadron neutral decay multiplicity: The charm tagging efficiency in the Monte Carlo is 
also observed to depend on the number of ir° mesons produced in the decay, even at fixed 
charged decay multiplicity, as the number of ir° mesons produced affects the amount of energy 
and transverse momentum available for the charged decay products. The average 7r° multiplicity 
in D + , D° and decays has been measured by MARK III p7j , but with large errors. However, 
for the D + (64±5) %, and for the D° (76.2±3.5) %, of the total decay width is to known exclusive 



final states, most of which have well measured branching fractions and low 7r° multiplicity [28] 

For the D+ and D°, the error is assessed by varying the Monte Carlo 7r° multiplicity in the 
unmeasured decay final states so as to reproduce the total variation in 7r° multiplicity allowed 
by the MARK III measurements, and adding a small contribution due to ir° multiplicity variation 
allowed by the measurement errors on the exclusive final state branching fractions. Since most 
of the decay modes with high tagging efficiency have well measured branching fractions, and the 
tagging efficiency for the unmeasured decay modes is lower than average, this procedure leads to 
total errors which are about half the size of those obtained by simply reweighting all decay final 
states to vary the 7r° multiplicity inclusively. The charged multiplicity and inclusive branching 
ratios to K° and A were held constant. The results of this procedure are systematic errors on 
R h of ±0.00016 for the D+ and ±0.00021 for the D°. For the D+, where only (18 ± 4) % of the 
total decay width is measured, the inclusive reweighting technique is used, leading to an error 
on i?b °f ±0.00016. The error due to charm baryons is negligible, due to their low production 
fraction and tagging efficiency. These errors are combined in quadrature to give a total error on 
R h of ±0.00030. 

Charm hadron to K° branching fraction: The inclusive branching ratios B(D + — > K°,K° ± X) 
B(B° -> K°,K0 + X), B(B+ -» K°,K° + X) and B(A+ -> A + X) were varied independently 
within the errors given by the Particle Data Group |26|], leaving the decay charged multiplicity 
distribution unaltered. The resulting variations in R^ were combined in quadrature to give a 
systematic error on i?b of ±0.00015. 

Charm semileptonic branching fraction: For semileptonic decays of charmed hadrons, an aver- 
age branching fraction of (9.3 ± 0.5) % was used. This value was obtained from two sources: the 
direct charm semileptonic branching fraction of (9.5 ±0.9) % measured at = 10 GeV [29|, and 



individual charm hadron branching fractions combined with the production fractions discussed 
above. The charm hadron semileptonic branching fractions were obtained from the branching 
fraction D° — ► e + X of (6.75±0.29) % |2q| , together with the measured charm hadron lifetimes |pq| 
and assumptions of lepton universality and equal semileptonic widths of all charm hadrons. This 
procedure gave a charm semileptonic branching fraction of (9.1 ± 0.6)%, which was averaged 
with the direct measurement. The direct measurement and D° — ► e + branching ratio uncertain- 
ties give a systematic error of ±0.00031 on R^ and an additional error of ±0.00015 results from 
the uncertainty in the charm hadron production fractions. 

Charm semileptonic decay modelling: The momentum spectra of the leptons in the rest frame 
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of the decaying charmed hadrons were modified according to the refined free-quark model of 



Altarelli et al. [3C]. The two parameters of the model, m s and pf, were chosen to be 0.001 GeV 
and 0.467 GeV, respectively, as given by a fit to DELCO @ and MARK III @ data per- 
formed by the LEP electroweak working group. Two sets of alternative values of the parameters, 
m s = 0.001 GeV, pf = 0.353 GeV and m s = 0.153 GeV, pf = 0.467 GeV, corresponding to the 
variation allowed by the fit, were used to estimate the systematic error on i?b of ±0.00029. 

Heavy quark production from gluon splitting: The production of heavy quark pairs via the 
processes g — > cc and g — > bb increases the tagging efficiency in charm and light quark 
events. The rate of g — > cc per multihadronic event has been measured by OPAL to be 
(2.38 ±0.48) x 10~ 2 Q, consistent with perturbative QCD calculations §. The rate of g — > bb 
has been measured by ALEPH and DELPHI j?4|] and has been averaged to give a value of 
(2.56 ±0.67) x 10~ 3 [ 2 1 1 , also consistent with perturbative QCD calculations. The Monte Carlo 



rates were adjusted to these central values, and the uncertainties in the rates lead to errors on 
R h of ±0.00018 from g -> cc and ±0.00027 from g bb. 

Inclusive K° and hyperon production: The total production rates of K°, A and other weakly 
decaying hyperons in the Monte Carlo were adjusted to agree with the values measured by 
OPAL |35| . The rates were varied by ±3.4%, ±6.5% and ±11.5% respectively, corresponding to 
the precision of the OPAL measurements combined with an additional uncertainty to take into 
account the extrapolation of the inclusive production rates to those for light quark events only. 
The results of these variations were combined in quadrature to give a systematic error on i?b 
of ±0.00001. The dependence of the folded vertex tagging efficiencies on the number of tracks 
produced in light quark events was found to be negligible. 

The finite number of charm and light quark Monte Carlo events contributes an additional error of 
±0.00010 on R h . 



8 Systematic errors: efficiency correlation 

The tagging efficiency correlation C b (as defined in section |2|) is determined from a sample of 7 million 
simulated bb events, with small corrections for the differing detector performance in each year. To 
evaluate the systematic error on this quantity, three classes of effect that can give rise to an efficiency 
correlation are considered: (1) kinematic correlations due to final state gluon radiation, (2) geometrical 
correlations due to detector non- uniformities, and (3) correlations coming from the determination of 
the primary vertex position. From equation ||| it can be seen that the fractional error on the correlation 
C b contributes directly to the fractional error on i?b, so an accurate evaluation of C b is essential. 

In general, correlations arise when the tagging efficiency e b depends on a variable or variables x, 
and the values of x are correlated between the two hemispheres of the event. The resulting tagging 
efficiency correlation C b can be evaluated as: 

b= (6 b (x)6 b (x)) 

(e h (x))(e\x)) 

where x and x are the values of x in the two hemispheres of the event, and the average is taken over 
all bb events in the sample. Such calculations, with the efficiencies evaluated in small bins of x, are 
used frequently in the correlation studies presented here. 

The total correlation estimates and errors are summarised in Table ||. The overall error on 
resulting from uncertainties in the correlation values is ±0.00066, corresponding to a relative uncer- 
tainty of 0.30 %. The systematic errors are evaluated by studying each component of the correlation 
separately. It is then checked that the total Monte Carlo correlation is reproduced by the sum of 
these components. The dependence of the Monte Carlo correlation on uncertainties in the detector 
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0.33 


0.22 
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0.03 


0.02 
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0.14 


0.12 




Data statistics 1993 


0.17 


0.11 




Data statistics 1994 


0.09 


0.06 




Data statistics 1995 


0.13 


0.09 


Primary vertex 


0.02 


0.02 


Simulation 


Detector resolution 


0.01 


0.02 




Beam spot size 


0.04 


0.06 




b quark fragmentation 


0.09 


0.08 




b hadron lifetime 


0.04 


0.03 




b decay multiplicity 


0.01 


0.01 


Total systematic error 


0.36 


0.25 


Monte Carlo statistics 


0.20 


0.17 


Total Correlation C b - 1 (%) 


1992 


0.84 ± 0.42 


1.01 ±0.33 




1993 


1.10 ±0.44 


1.13 ±0.33 




1994 


0.77 ±0.41 


0.92 ±0.31 




1995 


0.99 ± 0.42 


0.95 ±0.32 



Table 6: Errors on the hemisphere efficiency correlation C b — 1 for the vertex and combined tags, 
together with the correlation values and associated total errors used for each year of the data. The 
latter are obtained using the geometrical correlation values given in Table 0. 



simulation and physics modelling is also taken into account. This method differs from that used in 
where only relatively small Monte Carlo samples were available, and the overall correlation was de- 
termined by adding the estimates of each separate correlation component. Interdependence between 
the correlation components was treated in Q as an additional source of systematic error, but this is 
now accounted for by taking the correlation value directly from the Monte Carlo. 

The evaluation of the kinematic correlation is described in detail in Sections |3.1| an d |8T2| , and various 
comparisons between data and Monte Carlo are discussed in Sections [8J^ and |8.4|. The geometrical 
correlation, and the associated corrections to the Monte Carlo correlation value, are described in 
Section 8.5. Primary vertex effects are discussed in Section 3.6 and the detector simulation and 
physics modelling systematic errors in Section 8.7. The completeness test, that the total Monte 
Carlo correlation is consistent with the sum of the various components, is described in Section 3.5. 
Throughout this section, correlation values and uncertainties are given both for the vertex tag alone, 
and for the combined vertex and lepton tags. 



8.1 Kinematic Correlation — same-hemisphere events 

If an energetic gluon is radiated in a bb event, it may cause both b hadrons to recoil into the same 
thrust hemisphere. The tagging efficiency for such hemispheres is lower than that for the average b 
hemisphere, since the b hadrons have much lower momentum, and the double tagging efficiency for 
these events is much lower than for normal bb events since one hemisphere contains no b hadrons. 
These events therefore introduce a small hemisphere efficiency correlation, which can be calculated in 
the Monte Carlo from the number of such events and their hemisphere and double tagging efficiencies. 
The Monte Carlo predicts that 1.21 % of bb events passing the event selection have both b hadrons 
in the same hemisphere. However, such events constitute only 0.59 % of the tagged hemispheres and 
0.003 % of the double tagged events for the vertex tag, and only 0.62 % and 0.014 % for the combined 
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Figure 8: Distribution of the separation angle (f) vv between two vertices loosely tagged in the same 
hemisphere for genuine same-hemisphere bb events, normal opposite hemisphere bb events, charm and 
light quark events and data. 



tag. They contribute (0.02 ± 0.02) % to the tagging efficiency correlation for the vertex tag and 
(—0.03 ± 0.02) % for the combined tag, where the errors are due to Monte Carlo statistics. 

The number of same hemisphere events was compared in data and Monte Carlo by looking for 
events with two vertex tags in the same hemisphere, each passing a very loose neural network output 
cut of B > 0.5. In the Monte Carlo, 36% of these double tagged hemispheres come from genuine 
same- hemisphere events, 40 % from normal bb events with a mis-reconstructed vertex, and 23 % from 
charm and light quark events. The different contributions were statistically separated using the three- 
dimensional angle <p vv between the momentum vectors of the two vertices; the same-hemisphere bb 
events have a relatively broad 4> vv distribution and the remaining contributions are more concentrated 
at low values of <f> vv . The distribution of 4> vv for the double tagged hemispheres in data and Monte 
Carlo is shown in Figure ||. The 4> vv distribution in the data was fitted to the sum of the Monte Carlo 
distributions, allowing the rate of the Monte Carlo same hemisphere bb events to vary. The rate of 
genuine same hemisphere events was found to be consistent between data and Monte Carlo with a sta- 
tistical uncertainty of ±40 %. The systematic error due to uncertainties in the momentum distribution 
of the b hadrons was evaluated by comparing the predictions of JETSET 7.4 and HERWIG 5.9 [p^] . 
This gives a 7% relative uncertainty on the vertex tagging efficiency for same hemisphere events, 
and a 5 % uncertainty for the combined tag. The resulting total uncertainty on the size of the same 
hemisphere event correlation is ±0.08% for the vertex tag and ±0.06% for the combined tag. The 
tagging efficiencies are also sensitive to the modelling of b hadron decay in the Monte Carlo, but this 



is addressed in Section 8.7 



8.2 Kinematic Correlation — momentum and fragmentation 

After removing the same-hemisphere events in the Monte Carlo, each hemisphere of a bb event contains 
one b hadron. Kinematic variables, such as the momenta of the two b hadrons pb and pg, are correlated 
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between the two hemispheres due to final state gluon radiation. Since the b-tagging efficiency depends 
on these kinematic variables, such correlations can produce a tagging efficiency correlation. 

The b-tagging efficiency depends strongly on the b hadron momentum in the same hemisphere, 
since high-momentum b hadrons are likely to travel further before decaying, and decay producing 
higher momentum tracks which are measured with better resolution. Both of these effects produce a 
more easily resolvable secondary vertex. The efficiency is shown as a function of xb = PB/'-E'beam for 
the combined vertex and lepton tag in Figure ||(a). However, the b-tagging efficiency also depends on 
properties of the fragmentation tracks in the same hemisphere, for example their number, momentum 
and angular distribution, which influence the reconstruction of the primary and secondary vertices. 
The dependence of the b-tagging efficiency on the number of fragmentation tracks iVf rag is shown in 
Figure |9](b). Although -/Vf rag is negatively correlated withpB there is a dependence even at constant pb- 
As the number of fragmentation tracks in a hemisphere increases, the vertex tagging algorithm becomes 
increasingly likely to reconstruct the primary rather than the secondary vertex, and is therefore less 
likely to give a b-tag. 

Final state gluon radiation into one hemisphere decreases the momentum of the b hadron in that 
hemisphere, and increases the number of fragmentation tracks. It will also increase the hemisphere 



mass mh = y E^ 2 — p^ 2 where E^ and ph are the total energy and momentum in the hemisphere. Since 
the total momenta of the two hemispheres must balance, this will also reduce the total momentum, and 
hence the momentum available to the b hadron, in the opposite hemisphere of the event. Therefore 
the momentum pb of the b hadron in one hemisphere is correlated both with the momentum pg of 
the b hadron and the number of fragmentation tracks N-^g in the opposite hemisphere. This can 
be seen in Figure ^(c), where the Monte Carlo dependence of iVf rag in one hemisphere on pb in the 
other hemisphere is shown. At moderate values of pg, (iVf rag ) increases with increasing pg and it then 
decreases again at very high values of pg. 

As a result of these correlations, the b-tagging efficiency is found to depend not only on the same 
hemisphere b hadron momentum pb , but also weakly on the opposite hemisphere b hadron momentum 
pg , as shown in Figure ^(d) . This effect was also seen with the much simpler vertex tagging algorithm 
used in Q. The resulting kinematic correlation can be calculated by parameterising the b-tagging 
efficiency as a function of both pb and pg, and is given by: 

r b _ (e b (PB,P B ) eb (PBiPEi)) ... 
P ' P (<*(pB,Pt))(<*fa,PB))' U 

where e b (pi,p2) gives the b-tagging efficiency for a hemisphere containing a b hadron of momentum pi, 
the other hemisphere of the event containing a b hadron of momentum p2 , and the averages are taken 
over all bb events in the sample. The Monte Carlo predicts that the kinematic efficiency correlation 
is C£jj - 1 = (0.04 ± 0.05) % for the vertex tag and C* p - 1 = (0.06 ± 0.03) % for the combined tag. 
The size of the underlying momentum correlation between pb and pg is given by: 

r _ (pes Pb) 
PB ~ (pbXpb)' 

and the size of the resulting efficiency correlation, if the b-tagging efficiency depended only on the 
same hemisphere b hadron momentum, would be: 

r b = (e b (PB) e b (Pfi)) ^ 
P (e h (PB))(e h (p B )y U 

The Monte Carlo predicts a b hadron momentum correlation of Cp B — 1 = (0.80±0.01) %, in agreement 



with perturbative QCD calculations 37]. This would result in a tagging efficiency correlation of 



C p b - 1 = (0.87 ± 0.03) % for the vertex tag and C£ - 1 = (0.55 ± 0.03) % for the combined vertex and 
lepton tag. The additional dependence of the tagging efficiency on pg via the fragmentation tracks 
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Figure 9: Origins of kinematic correlations in Monte Carlo bb events: (a) Dependence of b-tagging 
efficiency e b on the scaled b hadron momentum xb = PB/'-E'bcam in the same hemisphere; (b) depen- 
dence of e b on the number of fragmentation tracks -/Vf rag in the same hemisphere, at various values of 
xb in the same hemisphere; (c) variation in the mean number of fragmentation tracks (-/Vf rag ) with xb 
in the opposite hemisphere; (d) dependence of e b on xb in the opposite hemisphere, at various values 
of xb in the same hemisphere. 
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reduces this correlation, to give an overall effect which is approximately zero. However, both the size 
of the momentum correlation Cp B and the separate correlation between -/Vf rag and pg are tested when 
calculating the systematic error. This is discussed in Section below. 

The momentum correlation has been treated above in terms of pb and iVf rag , quantities on which 
the tagging efficiency depends strongly, but whose hemisphere correlations are weak and difficult to 
measure. A different and complementary approach is provided by considering the total energy E^ and 
momentum p^ in each hemisphere, which are completely correlated by overall energy and momentum 
conservation in the event, and are affected by final state gluon radiation. These variables should 
give a good description of the kinematic correlation resulting purely from energy and momentum 
conservation. The variables E^ and ph are relatively easy to measure in the data, and the two 
hemispheres are related by E^ + E^ = 2 -Ebeam and ph = p^. However, the dependence of the tagging 
efficiency on these quantities is weaker than that on p B and -/Vf mg . 

In this approach, the kinematic correlation resulting from energy and momentum conservation 
may be calculated by parameterising the tagging efficiencies in terms of E^ and p^: 

b (e b (£ h ,p h )e b (£ h ,p h )) 

E « <e b (£ h ,p h )>(e b (£ h ,p £ )>' {) 

The Monte Carlo predicts that C% p - 1 = (-0.27 ± 0.04) % for the vertex tag and (-0.14 ± 0.04) % 
for the combined tag. If the efficiency is parameterised in terms of the hemisphere momentum ph 
alone, the correlation is (1.34 ± 0.03) % for the vertex tag and (0.85 ± 0.03) % for the combined tag, 
somewhat larger than the values obtained for C b using the b hadron momentum p B alone. 

The kinematic correlation C\ v — 1 is significantly more negative than Cp p — 1. To study whether 
significant additional correlation, other than that caused by energy and momentum conservation, is 
present in the Monte Carlo, the correlation Cp B between pb and pg was evaluated in bins of E^ and 
Ph- For the majority of the events, no significant additional correlation was observed, but for events 
with low values of ph, i.e. those with broad high mass jets in both hemispheres, a significant positive 
correlation between pb and pg was observed. This gives an extra contribution of about +0.15% to 
the tagging efficiency correlation, bringing the overall kinematic correlation closer to that calculated 
from pb and pg above. 

Both of the approaches described above were used to study correlations in the data. The correla- 
tions between pe and pg, and between -/Vf rag and pg, were checked using estimators of these quantities 



in a sample of loose double tagged events (see Section |8.3[) . The distributions of ph and E^, and the 
dependence of the tagging rates on them were studied. Monte Carlo events were reweighted to make 
the various distributions the same as those in the data, and the resulting changes in the correlation 
evaluated (see Section |8.4| ). Since all of the tests are somewhat indirect, and all the possible sources 
of kinematic correlations have not been completely understood, the largest systematic error resulting 
from any of the tests is used as the overall systematic error on the kinematic correlation. As can be 
seen from Table ^, the resulting systematic error is similar in size to the uncertainty in the correlation 
from Monte Carlo statistics. 

8.3 Kinematic Correlation tests — b momentum and fragmentation 

In order to study correlations in the data, a loose double tag was used to select a sample of almost pure 
bb events. Events were selected if each hemisphere contained a secondary vertex passing the neural 



network pre-selection described in Section 4.2 , i.e. at least three tracks with a decay length significance 
of at least \L/a\ > 3, or a lepton passing the selection described in Section [|. This selection has an 
efficiency for bb events of about 31 % and a b purity of about 95 %. 

Reconstructed secondary vertices were then used to form estimates p v of the corresponding b 
hadron momentum. The component of the b hadron's momentum associated with charged particles 
was estimated by summing the momenta of all charged particles in the jet, each weighted by X — the 
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weight for the track to have come from the b hadron decay described in Section 4.2. The neutral 
momentum component was estimated using the energy of all the electromagnetic calorimeter clusters 
in the jet, each weighted as a function of the angle between the cluster and the jet axis. If any charged 
tracks were associated to the cluster, their energies were first subtracted from the cluster energy, and 
the cluster was only used if the remaining energy was greater than zero. If the hemisphere was tagged 
only by a lepton, the lepton momentum was used as the p v estimator. The lepton momentum was 
scaled to give the same mean value of p v as that from reconstructed vertices. The correlation of 
x v = Pv/Ebeam with xb = PB/^beam is shown in Figure [K](a), and the distribution of x v in data 
and Monte Carlo loose double-tagged events is shown in Figure |lO|(b). The correlation coefficient 
P = ({pbPv) — (pb) (Pv)) / ( (J pb (J Pv) °f ^ e two quantities is 0.34, where ap B and <jp v are the standard 
deviations of pb and p v respectively. 

The number of tracks assigned to the hemisphere primary vertex -/Vp r i m was used as an estimator 
of the number of fragmentation tracks iVf rag in the hemisphere. The correlation of N pr [ m with iVf rag is 
shown in Figure |i~0|(c) , and the distribution of iV pr i m in data and Monte Carlo in Figure lC(d). The 
correlation coefficient of iVprim and iVf rag is 0.68. 

The correlation of the estimated b hadron momentum in one hemisphere, p v , with the estimated b 
hadron momentum in the other hemisphere, was measured to be Cp v — 1 = (0.670 ±0.023) % in the 
loose double tagged data sample, and Cp v — 1 = (0.664 ± 0.015) % in the corresponding Monte Carlo 
sample. The effects of including the same-hemisphere events and the non-bb background were checked 
in the Monte Carlo and found to be negligible. The Monte Carlo is therefore seen to reproduce the 
strength of the correlation between p v and p^ seen in the data. Since p v is only an estimate of pb 
the Monte Carlo was reweighted to change the true correlation between pb and pg and the resulting 
changes in Cp v and the efficiency correlation Cp studied. The statistical error on the difference between 
Cp v in data and Monte Carlo corresponds to an uncertainty of ±0.27% in Cp for the vertex tag and 
±0.19% for the combined tag. 

The size of the correlation between iVf rag and the opposite hemisphere b hadron momentum was 
also studied using the estimators described above. The average values of iV pr i m as a function of the 



estimated opposite hemisphere b momentum are shown in Figure 11(a) for the Monte Carlo and 



Figure |TT](b) for the data loose double tagged samples. A complicated correlation is seen in both data 
and Monte Carlo, with the value of (iV pr i m ) increasing in the region 0.2 < x v < 0.5 and then decreasing 
again in the region 0.5 < x v < 0.8. Qualitatively similar behaviour is also seen in the variation of 
(iVfrag) with xb shown in Figure ||(c). 

Two methods were used to compare the strength of the iV irag ,pg correlations seen in data and 
Monte Carlo. In the first method, the distributions were fitted to a parabola in the region 0.3 < x v < 
0.8 (where the bulk of the events lie), and the curvatures compared. The ratio between data and 
Monte Carlo curvatures was found to be 0.85 ± 0.10. In the second method, the root mean square 
deviation of (iV prmi )(x v ) from ((-/V prmi )) was calculated, weighted according to the distribution of x v . 
Here, (iV pr i m )(x v ) is the average value of iV pr i m at a particular value of x v , and ((iV pr i m )) is the average 
of iVprim calculated over all events. This method also takes into account the events outside the region 
0.3 < x v < 0.8, and gives a ratio between data and Monte Carlo of 0.91. The Monte Carlo was then 
reweighted to change the correlation between iVf rag and pg, and the variation of the difference Cpp—Cp 
studied. The full variation between the unweighted Monte Carlo sample and the sample reweighted 
to reduce the correlation between iV pr j m and p^ by 25 % was used to estimate the systematic error on 
Cpp — Cp, giving uncertainties of ±0.20 % and ±0.12 % for the vertex and combined tags respectively. 

The uncertainties on the separate components of the kinematic correlation Cp and Cpp — Cp were 
then added in quadrature to give total uncertainties on the kinematic correlation of 0.33 % and 0.22 % 
for the vertex and combined tags respectively. These uncertainties are the largest from any of the 
correlation tests, and set the size of the overall systematic error. The uncertainty for the vertex tag 
alone is somewhat larger for two reasons: the lepton tag has a weaker dependence on the b hadron 
momentum and no dependence on the number of fragmentation tracks; and the comparisons between 
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Figure 10: Estimators of of the scaled b hadron momentum xb and number of fragmentation tracks 
-^frag : (a) correlation of estimated xb (x v ) with true xb (b) distributions of x v in data (points) and 
Monte Carlo (histogram), (c) correlation of estimated Nf rag (iVp r im) and true iVf rag , (d) distributions 
of AWim hi data and Monte Carlo. 
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Figure 11: Estimated (A'frag) as a function of the estimated opposite hemisphere scaled momentum 
xb, in loose double tagged Monte Carlo (a) and data (b) samples. The fits used to estimate the size 
of the correlation are shown, (c) Hemisphere tagging efficiency correlations C b — 1 calculated using 
the estimated b hadron momentum function of p v only in the tagging hemisphere (circles), 

and in both the tagging and opposite hemispheres (squares). The data is shown by the filled symbols 
and the Monte Carlo by the open symbols. The errors are highly correlated between the points of 
different b-tagging efficiency, and the analysis cut is shown by the arrow. The individual points are 
slightly displaced along the x-axis for clarity. 
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data and Monte Carlo are dominated by the statistical precision of the tests, which are improved by 
adding the extra lepton tagged hemispheres. 

As an alternative approach, the distributions of p v and N pr i m were studied in data and Monte Carlo, 
and the influence of discrepancies in these distributions on the kinematic correlation investigated. In 
this approach, the decomposition of the kinematic correlation into the and Cpp — Cp components 
was not used. Event weights were calculated for the Monte Carlo to make the distributions of p v and 
iVprim match those in the data, and the resulting kinematic correlation calculated and compared with 
the value derived using: 

r b (e b (p v ,pv)e b (ffy,p v )) 

Pv ' m (e b (Pv,^))(e b (^,Pv)) 

i.e. the same as equation ||, but with p-Q and replaced by the estimators p v and p v and the averages 
taken over the sample of loose double tagged events. Event weights and efficiencies were calculated 
as a function of p v and p v , p v and iVp r i m , and -/Vp r i m and N^^, and the largest changes seen in C^ iP _ 
were 0.11 % and 0.07 % for the vertex and combined tags respectively, well within the systematic errors 
assigned for the kinematic correlation. 

Finally, the hemisphere tagging efficiency relative to the loose double tag pre-selection was calcu- 
lated directly in the data and Monte Carlo, from the fraction of events tagged as a function of p v , or 
of p v and p v . The two dimensional distribution of p v and p v can then be used to derive an effective 



tagging efficiency correlation. The results of this procedure are shown in Figure 11(c). The circles 
show the correlation calculated using the momentum dependence of the tagging efficiency in the tag 
hemisphere only: 



„ </(Pv)/(Pv)> 



c b (f( Pv )) 



</(Pv)></(p,)>' 

where f(p v ) is the fraction of hemispheres of estimated b hadron momentum p v that are tagged. The 
squares show the correlation calculated using the momentum dependence of the tagging efficiency in 
both hemispheres: 

ybftf- _ ^ _ </(Pv,l*)/(Pv,Pv)) 



C b (f( Pv ,Pv)) 



(/(Pv,Pv))(/(Pv,Pv))' 

In both cases, the correlations calculated from data and Monte Carlo are similar, and the increasing 
size of the correlation as the b-tagging efficiency is reduced is well modelled. The correlation using the 
tagging hemisphere momentum dependence alone is similar to that obtained in the Monte Carlo using 
the true b hadron momentum. The correlation due to the momentum dependence in both hemispheres 
is much closer to zero, again as in the Monte Carlo using the true b hadron momentum, though the 
statistical precision of this test is rather low. 

8.4 Kinematic Correlation tests — hemisphere energy and momentum 

To study the correlations expressed in terms of hemisphere energy and momentum, estimates were 
formed of the total energy and momentum in each hemisphere using both tracking and calorimeter 
information, and applying an algorithm to correct for double counting of charged particles [38]. These 



reconstructed energy and momentum distributions have widths of about 7GeV and are dominated 
by the experimental resolution. Since the hemisphere correlations of these quantities are known, the 
resolution of the estimators can be improved using information from both hemispheres and the beam 
energy constraint. This is done by calculating the two hemisphere masses and m£, using the beam 
energy to scale the total observed energy: 

earn / 7-1 2 9 

and then recalculating the hemisphere energies and momenta, using the two hemisphere masses, the 
known beam energy and two-body decay kinematics. The correlation of the reconstructed hemisphere 
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Tag 


1992 


1993 


1994 


1995 


Monte Carlo 


Vertex 


s~<h 

^geom 


" 1 (%) 


0.89 ±0.14 


1.15 ±0.17 


0.82 ±0.09 


1.04 ±0.13 


0.88 ± 0.02 


Lepton 


fib 

^geom 


" 1 (%) 


1.02 ±0.36 


0.54 ±0.30 


0.33 ±0.18 


0.12 ±0.29 


0.31 ±0.03 


Combined 


fib 

^geom 


" 1 (%) 


0.79 ±0.12 


0.91 ±0.11 


0.70 ± 0.06 


0.73 ± 0.09 


0.71 ±0.02 



Table 7: Geometrical hemisphere efficiency correlations Cg eom — 1 and their statistical errors, evaluated 
from each year of the data and from the Monte Carlo bb sample. 



mass with the true mass in Monte Carlo bb events is shown in Figure 12(a), and the reconstructed 
mass distributions for the loose double tagged data and Monte Carlo samples in Figure |i"2"Kb). 

The resulting distributions of hemisphere energies and momenta are shown in Figure [12|(c) and (d) . 
The resolutions of the estimators of hemisphere energy and momentum are 0.82 GeV and 0.68 GeV 
respectively, where the true quantities include unmeasured particles such as neutrinos. Some discrep- 
ancies are seen between the data and Monte Carlo distributions, which are addressed in the systematic 
error evaluation below. 

The systematic error on C\ v is evaluated in two steps. First, the two-dimensional distribution 
of hemisphere energies and momenta ph is reweighted in the loose double tagged Monte Carlo 
to match that of the data. This gives a shift in C\ v of (±0.08 ± 0.01)% for the vertex tag and 
(±0.04 ± 0.01) % for the combined tag, where the errors are due to data and Monte Carlo statistics. 
Secondly, the modelling of the tag efficiency as a function of E^ and p^ is checked, by studying the 
fraction of tagged hemispheres as a function of E^ and p^ in data and Monte Carlo, as shown in 



Figure 13. Reweighting the Monte Carlo to make the two dimensional distribution of the tag fraction 
as a function of and p^ agree results in a change in C\ of (±0.03 ± 0.15) % for the vertex tag and 
(±0.03 ± 0.13) % for the combined tag, where the errors are statistical. These changes are all small, 
and well within the overall systematic errors assigned for the kinematic correlation. 



8.5 Geometrical Correlation 

The two b hadrons in a bb event tend to be produced back-to-back, so their decay products are likely 
to hit geometrically opposite parts of the detector. This introduces an efficiency correlation if the 
tagging efficiency is not directionally uniform. The vertex tag efficiency depends on the polar angle 
of the b hadron because multiple scattering degrades the tracking resolution more as cos 6 increases. 
For some of the data, the tagging efficiency also depends on </>, due to inefficient regions of the silicon 
detector. The lepton tag efficiency depends on both and (f>, due to localised regions of inefficiency 
caused by effects such as incomplete muon chamber coverage. 

The size of the correlation from this effect can be estimated from the data, by measuring the 
hemisphere tagging probability as a function of the thrust axis polar angle 6t and azimuthal angle 

, b . Mf + (0r,<h)f-(0T,<h)) 



a 



gcom 



(f+(0 T , 0r) + f-(9 T , cf> T ))(f+(6 T , cfr) + f-(9 T , <f> T )) (?) 

where / + (/ _ ) is the fraction of hemispheres in the +z (—z) direction that are tagged, and the averages 
are performed over all values of cos 9t and 4>t- The values of / + and /~ are calculated in small bins 
of cos 0t and (f>T- As the tagged samples are very pure in bb events, contributions from charm and 
light flavour events to the correlations calculated from the tagged data are negligible. 

The resulting geometrical correlations for the vertex, lepton and combined tags are given for each 
year of the data in Table ^. The values for the different years are consistent, and in agreement with 
those derived from Monte Carlo simulation of the detector configuration in each year. 

In order to calculate the overall correlation coefficient appropriate for each year of the data, the 
geometrical correlation in the Monte Carlo sample was calculated using equation [5]. These values are 
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Figure 12: (a) Correlation of reconstructed and true hemisphere masses m^; (b-d) Distributions of 
reconstructed hemisphere mass, energy and momentum in loose double tagged events in data (points) 
and Monte Carlo (histogram). 
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Figure 13: Fractions of tagged hemispheres in loose double tagged events as functions of (a) hemisphere 
energy and (b) hemisphere momentum p^. The data are shown by the points and Monte Carlo by 
the histograms. 



also given in Table 0. A correction: 



AC 1 



/"<b,data 
^geom 



c. 



b,MC 



(8) 



was then applied to the overall Monte Carlo correlation to give the final correlation value for each 
year of the data shown in Table [6|. 

The geometrical correlation is calculated using the thrust axis as an estimator of the b hadron 
direction, assuming the two b hadrons to be back to back. Final state gluon radiation can push the 
two b hadrons away from the thrust axis direction, and towards regions of the detector with different 
b tagging efficiency. The size of the true geometrical correlation is therefore given not by equation |7| 
but by: 

(e b (cos 6> B , <Pb ) e b (cos 6>g , 4> B ) ) 



ct, 



(9) 



ge ° m <6 b (cOS#B,0B)}(e b (cOS%,0B)>' 

where cos#B) <t>B and cos#g, </>g give the directions of the two b hadrons, approximated by the thrust 
axis direction in equation [7|. The error due to this approximation was estimated by calculating AC b 
using equation |8[ but with replaced by C^j^p calculated from an independent Monte Carlo 

sample with a geometrical correlation twice as strong as that seen in the data samples. The value of 
AC b calculated using the thrust axis direction was then compared with that calculated using the true 
b-hadron directions, evaluating the two C b <^p values using equation [9|. The difference was 0.03% for 
the vertex tag alone and 0.02 % for the combined tag, which was taken as the systematic error on the 
geometrical correlation correction. 

8.6 Primary Vertex Correlation 

In the previous analysis Q, the primary vertex position was determined using all tracks in the event, 
and then used in the tagging of each hemisphere. This led to a small but significant negative efficiency 
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correlation. In this analysis, the primary vertex correlation has been largely eliminated by determining 



the primary vertex separately for each hemisphere (see Section |4.l| ). However, the beam spot constraint 
is still shared between the hemispheres, which may lead to a small residual correlation. 

This was investigated in the Monte Carlo by generating separate independent beam spot constraints 
for each hemisphere, according to a Gaussian distribution of the same width as the usual beam spot 
constraint, centred about the true event primary vertex position. This procedure retains the constraint 
from the beam spot but removes any correlation due to the shared position information. The Monte 
Carlo tagging efficiency correlation using this modified constraint was found to be different by only 
0.02 % from that with the normal constraint, showing there is no significant primary vertex correlation 
from the shared beamspot. The full size of this shift is taken as an additional systematic error on the 
overall correlation value. 

8.7 Correlation Simulation Systematics 

The estimate of the overall size of the efficiency correlation from the Monte Carlo relies on the detector 
simulation and the modelling of the production and decay of b hadrons. The following uncertainties 
were considered to estimate the systematic error from these sources. The effect of each variation on 
the correlation values is given in Table ^. 



Detector Resolution: The detector resolution was varied using the method discussed in section |7.1| , 
varying the r-<j) resolution parameter in the range 0.9-1.1. Since the effect seen is very small, 
additional studies varying the r-z tracking resolution, tracking efficiency, silicon alignment and 
silicon hit matching efficiency were not carried out. 

Beam spot size: The size of the LEP beam spot varied between each year of the data due to the 
differing operating conditions of the accelerator. The beam spot size was measured in each year 
using Z° — > events, and was typically 140 fim in x, 10 /xm in y, and 7 mm in z. In x and z, 

the beam spot size also decreases gradually during the course of each fill, due to the operation 
of the LEP wiggler magnets. Residual uncertainties in the beam spot size are estimated to 
be ±10 /um in x, ±5/um in y and Umm in z. The effective beam spot constraint in Monte 
Carlo events was varied by these amounts to assess the resulting uncertainty in the correlation 
coefficients. 

b quark fragmentation: The b quark fragmentation was varied in the Monte Carlo by applying a 



weight to each simulated event using the fragmentation functions discussed in Section |7.4| , so 
as to vary the average scaled energy {xe} of the weakly decaying b hadrons by ±0.008. This 
variation represents the accuracy of {xe} measured by LEP experiments 



b hadron lifetime: The lifetimes of the b hadrons were varied simultaneously by ±0.05 ps using a 
weighting method. The size of the variation was chosen to be larger than the accuracy of the 
world average b-lifetime |2(| to allow for the uncertainties due to different efficiencies for different 
b hadron species. 

b hadron charged decay multiplicity: The average charged decay multiplicity of the b hadrons 
was varied by ±0.35 using a weighting method. The size of the variation reflects the accuracy 
of the measurements by DELPHI and OPAL 



The resulting variations in the correlation coefficients are all rather small, and much smaller than the 
dominant uncertainties caused by momentum correlation effects. 



36 



ftl.75 

~ 1.5 

U 

cl.25 
© 

'X 1 

gO.75 
3 0.5 
0.25 




~i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — | — i — i — i — 
□ Geometrical • Overall correlation fxn \ J 

A Same hemisphere O Component sum vJx 
Momentum 



no 



rf> 



A* A° A° A° 



A v 



A 



,0 



.0 



A v A A v 
_i i I i i i I i i i_ 



18 



20 



22 



24 



26 28 30 

b-tag efficiency (%) 



Figure 14: Comparison of the true overall tagging efficiency correlation and the various components 
for the combined vertex and lepton tag, as the b-tagging efficiency is varied by changing the vertex 
tag cut. The individual components are slightly displaced along the x-axis for clarity, and the analysis 
cut is shown by the arrow. 
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Vertex 


Combined 


Same hemisphere events 


0.02 
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0.02 


-0.03 


±0.02 


Momentum correlation 


0.04 
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0.05 


0.06 


±0.03 


Geometrical correlation 


0.88 
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0.02 


0.71 


±0.02 


Component sum 


0.94 


± 


0.06 


0.74 


±0.04 


Overall correlation 


0.83 


± 


0.20 


0.93 


±0.17 



Table 8: Components of the hemisphere efficiency correlation C b — 1 in bb Monte Carlo events, 
together with their sum and the overall correlation in the Monte Carlo sample. Only statistical errors 
are included. 



8.8 Correlation Completeness Test 

The sizes of the various components of the correlation in Monte Carlo, together with their associated 
statistical errors, are given in Table |8[ The evolution of the various components as the b-tagging 
efficiency is varied is shown in Figure [l4|. In the absence of interdependence between the correlation 
components, the values of C h — 1 from each component can be added to give the overall correlation. The 
sum of the components in Table |8| is compatible with the overall correlation evaluated from the ratio of 
the event tagging efficiency to the product of the hemisphere tagging efficiencies, for both vertex and 
combined tags. This strongly suggests that no large correlation contribution has been omitted from 
the systematic error evaluation, though the precision of this test is limited by the overall Monte Carlo 
statistics and by interdependence between the different components. In the previous analysis H, an 
additional systematic error was included to cover such interdependence, since the overall correlation 
was determined by simply adding the components. No such error is included here, since the Monte 
Carlo is used to determine the value of the correlation directly, including any interdependence. 
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These correlation component sums have been made using the C^p estimate from equation || for the 
momentum correlation. If the energy and momentum conservation C\ estimate from equation ^ is 
used instead, the component sums are (0.63 ±0.05) % and (0.54 ±0.04) % for the vertex and combined 
tags, still compatible with the overall Monte Carlo correlation. It should be emphasised that the central 
value used for the correlation does not depend on which method is used to evaluate the momentum 
correlation component. 



9 Summary and Conclusion 

The fraction of Z° — > bb events in hadronic Z° decays was measured using data collected by OPAL 
between 1992 and 1995 using tags based on displaced secondary vertices and identified leptons, giving 
the result 

R h = 0.2178 ± 0.0011 ± 0.0013 

where the first error is statistical and the second systematic. The systematic error does not include 
the effects of varying R c from its Standard Model expectation. The result depends on R c as follows: 

AR h AR C 
-IT = -0.056—^ 
/lb K c 

where AR C is the deviation of R c from the value 0.172 predicted by the Standard Model. 

This result is in agreement with and supersedes our previous measurement [||. It is also in agree- 
ment with other recent measurements at LEP ]l2|, p|] and SLC [42], and with the Standard Model 
prediction of 0.2155 ± 0.0003 @. 
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